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Abstract: Correlation functions of operators with a conformal dimension of 0(N 2 ) 
are not well approximated by the planar limit. The non-planar diagrams, which in the 
bulk spacetime correspond to string loop corrections, are enhanced by huge combinato- 
rial factors. In this article we show how these loop corrections can be resummed. As a 
typical example of our results, in the half-BPS background of M maximal giant gravi- 
tons we find the usual 1/N expansion is replaced by a 1/ (M + N) expansion. Further, 
we find that there is a simple exact relationship between amplitudes computed in the 
trivial background and amplitudes computed in the background of M maximal giant 
gravitons. Finally, we also find strong evidence for the BMN-type sectors suggested in 
arXiv:0801.4457. The decoupling limit of arXiv:0801.4457 captures the decoupled low 
energy world volume theory of the intersecting giant graviton system and this theory 
is weakly coupled even when the original M = 4 super Yang-Mills theory is strongly 
coupled. 
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1. Introduction 

According to the AdS/CFT correspondence [1], the conformal dimension of an operator 
in the J\f = 4 super Yang-Mills theory maps into the energy of the corresponding state 
in IIB string theory on the AdS 5 xS 5 background. Thus, operators with a very large 
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dimension will map into states with a very large energy. If this energy is large enough, 
back reaction can not be neglected and the state is best thought of as a new geometry 
which is only asymptotically AdS 5 xS 5 . Good examples of such operators include the 
Schur polynomials [2, 3, 4] with 7£-charge of 0(N 2 ) (dual to LLM geometries[5, 6]) and 
the operators obtained by distributing a gas of defects on Schur polynomials (which 
seem to be dual to asymptotically AdS 5 xS 5 charged black holes[7, 8]). To build a 
proper understanding of these operators in the gauge theory one would like, at least, 
to compute the anomalous dimension of these operators and to deal with their mixing. 

In this article we consider the problem of computing correlators of heavy operators. 
The operators we have in mind are a small perturbation of a BPS operator which 
has 7£-charge ~ conformal dimension ~ N 2 . In the language of [9] we study almost 
BPS operators. To solve this problem one must reorganize perturbation theory by 
resumming an infinite number of diagrams. The need for this reorganization is that 
for these operators non-planar diagrams can't be neglected [10]. A similar situation 
is provided by the BMN sector [11] of M = 4 super Yang-Mills theory. In this case 
one considers operators with 7£-charge J ~ 0(VN). One finds the usual expansion 
parameter is replaced by a new expansion parameter equal to ^[12]. Thus, one must 
take <C 1 to suppress non-planar diagrams. 

The question we ask and answer in this paper is: 

Is there a reorganization of the expansion for almost BPS operators which have 
IZ-charge of 0(N 2 ) and if so, what is the expansion parameter? 

We focus on almost BPS operators since we want to extrapolate our computations to 
strong coupling, where we can compare with the dual string theory. 

A particularly useful way to describe the half-BPS sector is to use the Schur 
polynomials [2, 3, 4]. Among the operators we consider, are small perturbations of 
a Schur polynomial Xb(Z) with B a Young diagram that has M columns and N rows, 
and M is 0(N). This corresponds to an LLM geometry [5] with boundary condition 
that is an annulus. The inner radius of the annulus is oc \J~M and the outer radius is 
oc \J M + N. We demonstrate, in section 2.1, that there is a reorganization of the per- 
turbation theory in the half-BPS sector and that the new expansion parameter is M + N - 
In section 2.2 we confirm this answer by studying the holography of the IIB supergrav- 
ity in the relevant LLM geometry. We then generalize these results to multi-ring LLM 
geometries (in section 2.3), backgrounds with more than one charge (in sections 2.4 and 
2.5) and beyond the BPS sector (in section 3). The Schwinger- Dyson equations provide 
a very powerful approach to the computation of correlators in the annulus background. 
We present these details in the Appendices. 
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2. Half-BPS Sector 



Schur polynomials provide a very convenient reorganization of the half-BPS sector. 
This is due to the fact that their two point function is known to all orders in -^[2, 3] 
and that they satisfy a product rule allowing computation of exact n-point correlators 
using only two-point functions (see Appendix D for a summary of the results we use). 
In this section we will use this Schur technology to provide an answer, in the half-BPS 
sector, to the question posed above. 

M = 4 super Yang-Mills theory has 6 scalars fa transforming in the adjoint of 
the gauge group and in the 6 of the SU{A)n symmetry. We shall use the complex 
combinations 

Z = 0i + ifa, Y = 3 + 204, X = 05 + 206, 

in what follows. We focus on the contributions coming from the color combinatorics; 
we drop all spacetime dependence from two point correlators 1 , that is, we use the two 
point functions 

(Z y 4) = (y^) = (XtjXl) = 5aS jk . (2.1) 

Since we are not explicitly displaying the spacetime dependences, it is important to 
point out that all holomorphic operators are inserted at a specific spacetime event and 
all anti- holomorphic correlators are inserted at a second spacetime event. These correla- 
tors are called extremal correlators; there are non-renormalization theorems protecting 
these correlators [13, 14]. 

When we talk about a half-BPS operator in what follows, we mean an operator 
built only from Zs. These operators will not break any further supersymmetries beyond 
those broken by the background itself. We obtain almost BPS operators by sprinkling 
Ys and Xs in the operator. 

2.1 Super Yang-Mills Amplitudes 

The exact computation of multi-trace correlators at zero coupling is most easily carried 
out by expressing the multi-trace operators of interest in terms of Schur polynomials 

J]Tr (Z*) = $> flX fl(Z), n Tr (( Zt ) mj ) = E/W^)- ( 2 - 2 ) 

i R j R 

The coefficients and /3r appearing in the above expansion have no dependence on 
N. It will be useful first to compute these correlators with a trivial background. A 

*It is simple to reinstate the spacetime dependence in the final result. 
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little Schur magic now gives (see Appendix D) 

Adm^N) = ll[Tr(Z^)Tr((Z^) \ = £ a R p T ( Xr {Z) Xt {Z^)) 

\ i,j I R,T 

= $> R /^- (2.3) 

R 

In this last equation f R is the product of the weights 2 (one for each box in the Young 
diagram) for Young diagram R. 

In what follows, B denotes the rectangular Young diagram with N rows and M 
columns. We will often refer to B as the annulus background because the corresponding 
LLM geometry is obtained by taking the LLM boundary condition to be a black annulus 
on a white plane. The expectation value of an operator O in background B is given by 

_ (xb(Z)xb(Z')Q) 
{U) "- (xb(Z)xb(Z>)) ■ (2A) 

Our normalization is chosen so that the expectation value of the identity is 1. Recall 
that B is a Young diagram with N rows and M columns. The product rule satisfied 
by the Schur polynomials can be viewed as a consequence of the fact that the Schur 
polynomials themselves are characters of SU(N) and that (i) the character of a direct 
product of representations is just the product of the characters and (ii) the character 
of a given reducible representation is equal to the sum of characters of the irreducible 
representations appearing in the reducible representation. In the product 

Xr(Z)xs(Z) = ^29rstXt(Z), 

T 

the Littlewood-Richardson number g^sr counts the number of times irreducible rep- 
resentation T appears in the direct product of the irreducible representations R and 
S. Since we choose our background B to be a Young diagram with M columns and N 
rows, it is a singlet of SU(N) and consequently we have the product shown in figure 1. 
Using this product rule, the correlator in the annulus background is 

A B ({n, m ,},N) = (nTr(Z»')T r ((Zt r A = ^ aA (xMx-Mx*{*)Xt&)) 

\ i,j I B R > T 

\- R ( x+ r(z)xM&)) 

= l^ a ^T J- 

R,T JB 

R iB 



2 Recall that the box in row i and column j has a weight N + j ■ 
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B 



x R 



= +R 



Figure 1: The product rule used to compute correlators in background B. This figure defines 
+R. 

Recall that f +R is the product of the weights appearing in Young diagram +R and 
J'b is the product of the weights appearing in Young diagram B. All of the weights 
appearing in the product f B are repeated in the product f +R so that after canceling 
common factors J+r/ f B is simply equal to the product of the weights of the extra boxes 
stacked to the right of B to form +R. Thus, /+#/ f B is obtained from f R by replacing 
N —> N + M . Comparing (2.5) and (2.3) we find 

A B ({rn; rrij}, N) = A({ni, rrij}, M + N). (2.6) 

We know that the correlator A({rii, rrij}, N) admits an expansion in (2.6) tells us 
that A B ({rii] raj}, N) admits an expansion in N + M - If we assume that Yl,i n i ~ 0(1), 
we obtain correlators of operators that are dual to gravitons. Thus, for gravitons in 
the background B N ^ M clearly plays the role of a loop expansion parameter. Thus, 
we learn that there is a reorganization of the jj expansion for these correlators in the 
annulus background B and further, that the new expansion parameter is N } M - Our 
relation (2.6) implies that the only effect of the background is to shift N — > N + M. 

One can easily check that this relation (2.6) is not a property of the full theory. 
Although a slight modification of (2.6) does allow us to relate the one point functions 3 
(Tr (Z n Z J<n )) and (Tr (Z n Z^ n )) B , we have not worked out a relation between amplitudes 
in general. Our interest in the one point functions (Tr (Z n Z^ n )) B is that they appear 
in the intermediate steps of the computation of correlators of BMN-like probes of the 
annulus so that we manage to obtain a simple relation between the trivial background 
and the annulus background for both the BPS and near-BPS sectors of the theory. 

The amplitudes A({ni] rrij},M + N) are the amplitudes of a theory with a gauge 
group of rank N+M and no background. In the LLM language, the boundary condition 

3 Again computed in the theory with gauge group U(N + M). 
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for this geometry is simply a black disk of radius yN + M . In our theory (with gauge 
group of rank N) the background B corresponds to an annulus with inner radius \f~M 
and outer radius \JM + N. Thus, one way to interpret the relation (2.6) is that any of 
the half-BPS probes considered above are unable to detect the hole in the middle of the 
annulus. We have not considered probes built using J? instead of Z; these will detect 
the hole. Indeed, acting with traces of Z on the background Xb{Z) produces a new 
Schur polynomial that has extra boxes stacked adjacent to the upper right hand corner 
of B - this corner maps to the outer edge of the annulus defining the LLM boundary 
condition. Acting with erodes corners from the lower right corner of B - this corner 
maps to the inner edge of the annulus. See [16] for further details. Probes built using 
22 instead of Z are also half-BPS probes. 

We have related amplitudes in the gauge theory with a background of M giant 
gravitons and gauge group U (N) to amplitudes in the gauge theory with trivial back- 
ground and gauge group U (N+M). This is reminiscent of the infrared duality proposed 
in [15] which exchanges the rank of the gauge group and the number of giant gravitons. 
In fact, (2.6) tells us that half-BPS correlators computed in the U(N) gauge theory in 
the background of M giant gravitons are exactly equal to the same half-BPS correlators 
computed in the U(M) gauge theory in the background of N giant gravitons. Since 
these correlators are extremal and hence not renormalized, our computations give the 
value of these correlators in the deep infrared limit of the gauge theory and thus seem 
to provide nontrivial support for the duality proposed in [15]. Note however, that cor- 
relators of operators built using J? (which are also extremal) will not agree. This is 
not obviously in conflict with the proposed duality of [15]. In [15] near extremal black 
holes are considered. These can be understood as a condensate of giant gravitons [17] 
represented in M = 4 super Yang-Mills theory by a Schur polynomial corresponding 
to a triangular Young diagram [6]. The LLM boundary condition for our background 
is a black annulus and ^ correlators explore the inner edge of the annulus. The LLM 
boundary condition dual to the Schur polynomial with triangular Young diagram label 
is a gray disk; there is no inner edge. It would be interesting to see how many of our 
results for the black annulus boundary condition can be generalized to the gray disk 
boundary condition. This generalization is nontrivial. 

The main evidence given in [15] for the proposed duality was an exact correspon- 
dence between the gravitational entropy formulae of small and large charge solutions. 
We can give arguments, in the free field theory, that suggest that the entropy of the 
state formed by N g condensed giant gravitons in the theory obtained by taking the 
near-horizon limit of N D3-branes is equal to the entropy of the state formed by N 
condensed giant gravitons in the theory obtained by taking the near-horizon limit of N g 
D3-branes. This provides further evidence for the duality of [15] in a completely differ- 
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ent region of parameter space to that probed by supergravity. According to [7] the near 
extremal black holes in AdS space can be obtained by attaching open string excitations 
to the Schur polynomials with triangular Young diagram labels. If we are very close to 
extremality, the number of open string excitations attached is much smaller compared 
to the total number of fields in the operator. In this situation instead of attaching open 
strings by adding boxes to the original triangular label, it should be a good approx- 
imation to fix the tableau shape and replace some boxes (in arbitrary places on the 
Young diagram) with open strings. Assuming that the state obtained by deleting the 
open strings is again a typical state, i.e. again a triangle we have something like the 
situation given in figure 2. To compute the entropy associated with these states, we 
need to count the number of different ways of attaching the open string defects. This is 
given by the number of different ways we can pull the occupied boxes off the triangular 
Young diagram. We don't know exactly how to compute this number. However, it is 
easy to argue that this number is invariant under flipping the Young diagram so that 
the two states in figure 2 are swapped: this follows from two facts (i) the counting 
problem is constrained by the rule that we can pull boxes off in any order as long as 
after each box is pulled off we continue to have a valid Young diagram 4 and (ii) a shape 
that is (is not) a valid Young diagram before the flip is (is not) a valid Young diagram 
after the flip. Thus, the entropy of these two states agree. 



Figure 2: The typical state of N g condensed giant gravitons in the theory with gauge group 
U (N) (on the left) and the typical state of N condensed giant gravitons in the theory with 
gauge group U(N g ) (on the right). The black stripe represents the boxes occupied by open 
string defects. 

A few comments about the above flip of the Young diagram are in order. The 
half-BPS sector of M = 4 super Yang-Mills theory can be mapped to the quantum Hall 
system with filling factor equal to one [22, 23]. In the half-BPS sector the above flip 

4 This is required because removing each box must define a valid subduction. See[18, 19, 20, 21]. 
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of the Young diagram is a Z 2 symmetry that exchanges particles and holes[22, 23, 24]. 
If one employs a continuum (field theory) description of the quantum Hall fluid, the 
resulting fluid is incompressible so that the corresponding continuum Lagrangian has 
a gauge invariance under area preserving diffeomorphisms[25]. Small fluctuations of 
the density, for a fluid of charged particles in a background magnetic field, are well 
described by a Chern-Simons theory; the U(l) gauge invariance is nothing but the 
area preserving diffeomorphisms for the small fluctuations [25]. The fluid description 
correctly captures the long distance physics of the quantum Hall effect. It does not 
capture the fact that, since the fluid is described by N electrons, it has an intrinsic 
granular structure. One can capture this granular structure by using a Chern-Simons 
matrix model description [25]. In this Chern-Simons matrix model description the above 
Z 2 symmetry is nothing but level/rank duality of the Chern-Simons matrix model[26, 



One puzzling feature of the duality of [15] is the fact that it mixes giant gravitons 
and the branes (called background branes in [15]) whose near horizon geometry is the 
AdS space we work in. There is a big difference between the background branes and the 
giant gravitons. Indeed, in string theory the background branes carry a net RR-charge; 
giant gravitons carry no net charge - they are dipoles. In the dual super Yang-Mills 
theory, changing the number of background branes changes the rank of the gauge group 
i.e. the number of fields we integrate over when performing a path integral quantization. 
Changing the number of giant gravitons leaves the rank of the gauge group unchanged, 
but it does change the background i.e. the integrand we use when performing a path 
integral quantization. If the duality of [15] is correct, we need to understand how, in 
this case, changing the integrand has exactly the same effect as changing the number 
of variables over which we integrate. Computing |-BPS correlators in the annulus 
geometry gives a nice toy model in which to explore this issue. This is because we can 
reduce the whole problem to eigenvalue dynamics in zero dimensions. Consider the 
computation of the correlator (Tr (Z)Tr (^)) in the trivial vacuum of the U(N + M) 
theory (A is usual the Van der Monde determinant) 



The same correlator in the U (N) theory with a background of M giant gravitons is 



We already know that (2.7) and (2.8) give the same result. The reason why the two 
agree is now evident: in (2.7) we integrate over an extra M variables; these extra 



27]. 




(2.7) 




(2.8) 
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contributions add to give a larger result than that obtained for the same correlator in 
the U(N) theory. In (2.8) the extra factor in the integrand implies that the integrand 
now peaks at larger values for the \zi\; this again gives a larger result than that obtained 
for the same correlator in the U(N) theory in the trivial vacuum. Our computation 
gives a simple picture of how, in this case, changing the integrand has exactly the same 
effect as changing the number of variables over which we integrate. It also suggests 
that the duality of [15] might be derived by starting with a suitable U (N + N g ) theory 
and (i) integrating out the degrees of freedom associated with N g colors to get a U (N) 
gauge theory with a background given by a condensate of N g giant gravitons or (ii) 
integrating out the degrees of freedom associated with N colors to get a U(N g ) gauge 
theory with a background given by a condensate of N giant gravitons. 

2.2 Supergravity Amplitudes 

According to the AdS/CFT correspondence, correlation functions can be computed in 
the strong coupling limit of M = 4 super Yang-Mills theory using the dual hologram, 
which is type IIB supergravity. The paper [28] , has given a powerful general approach to 
holography in the LLM backgrounds, generalizing and extending the Coulomb branch 
analysis of [29, 30]. The formalism given in [28] is an application of the general method 
of [31] which employs the method of holographic renormalization[32]. An important 
result of [28] was the demonstration that the asymptotics of the LLM solutions correctly 
encode the vacuum expectation values of all single trace | BPS operators to the leading 
order in the large N expansion. Can we reproduce the results of the last section using 
holography in the LLM background? 

One case of interest to us in this article is that of graviton three point functions, 
in the background created by a heavy operator. We have computed these correlation 
functions in the free field theory; supergravity will reproduce these correlators in the 
strong coupling limit. We will now argue that it is natural to expect that these two 
results will agree since we can argue that there is a non-renormalization theorem pro- 
tecting the three point functions of interest to us. The logic of this argument is very 
similar to the argument of [28] which argued that the generic one point function in the 
LLM background is protected. The three point functions we are interested in 

(Tr(Z n )Tr(Z m )Tr(Z tm+n )) B 

can of course, also be understood as a ratio of (higher point) correlators in the original 
trivial background 

( XB (Z) X B(Zr)Tr(Z n )Tr(Z m )Tr(& m+n )) 
( XB (Z) XB (&)) 
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The correlators appearing in the above expressions are extremal correlators in the 
language of [33]. The computations of [33] show that at the leading order in large 
N, the extremal correlators take the same value at large 't Hooft coupling as in the 
free field theory and hence it is natural to expect that all extremal correlators are not 
renormalized 5 . 

To extract three point functions using the methods of [28], we would need to solve 
the equations of motion, to quadratic order, around the LLM solution. We have not 
managed to do this. For the case of extremal correlators, there is some room for 
optimism: the analysis of [33] argued that for extremal correlators, the bulk cubic 
supergravity coupling vanishes and the entire contribution comes from a boundary 
term. The LLM geometries are asymptotically AdS so that one might have hoped that 
there was a simple explanation of the results of the previous section. We have not found 
one. Note also that the fact that the bulk cubic supergravity coupling vanishes in the 
AdSs background need not imply that it continues to vanish in the LLM background. 

Rather than pursuing the computation of three (and higher) point functions di- 
rectly, we have found it simpler to relate the three point functions we'd like to compute 
to one point functions, since these have already been obtained in [28]. For concreteness, 
consider the three point function (B is again the annulus diagram) 

(Tr(Z 2 )Tr(Z 2 )Tr(Z t4 )) B 

which can also be written as the one point function of Tr(Z 2 ) in the normalized state 

|$)=AT(Tr(Z 2 ) + Tr(Z 4 )) XB (Z)|0). 

Af is a normalization factor. Now, summing only planar diagrams 

(Tr(Z 2 )Tr(Z 2 )Tr(Z t4 )) = 16N 3 

so that according to our result of the previous section 

(Tr(Z 2 )Tr(Z 2 )Tr(Z t4 )) B = 16(N + Mf . 

The supergravity one point function is computed with a normalized state |<E>); the above 
three point function is computed with a different normalization. Noting that (we are 
again using results from the previous section) 6 

((Tr (Z 2 ) + Tr (Z 4 ) )(Tr {Z^) + Tr (Z^)) Xb (Z) Xb (Z^)) 

5 Note that the supergravity result suggests that the planar contribution is not renormalized. Here 
we are using the stronger conjecture [33] which claims the non-renormalization for any N. Our results 
provide further evidence for this stronger conjecture. 

6 Note that since |$) has two additive components the one-point function has additional cross terms; 
these however evaluate to zero leaving only the term we wish to compute. 
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= [4(N + Mf + 22(iV + M) 2 ] f B = 4(N + M) 4 f B (l + 0((N + My 2 )) , 
we easily find 

($|Tr(Z 2 )|$> = — -i— . 
v 71 ' N + M 

This is the result we would like to reproduce from supergravity. 

We will summarize the main steps involved in extracting correlators of local oper- 
ators in the boundary gauge theory, from the ten dimensional asymptotically AdSsxS 5 
solutions of Type IIB supergravity. For all the details and further references consult 
[28]. Given an asymptotically AdSsxS 5 solution of Type IIB supergravity, one must 
first systematically reduce the ten dimensional solution to a solution of five dimensional 
gravity coupled to an infinite number of five dimensional fields. If this reduction is per- 
formed in terms of gauge invariant variables (as it is in [28]), it does not matter in which 
gauge the original solution is supplied. This reduction is then supplemented with a non- 
linear field redefinition, performed in such a way that the equations of motion in terms 
of the redefined fields can be obtained by minimizing a local five dimensional action. 
One can then apply the usual holographic rules to compute the correlation function of 
boundary operators from the bulk five dimensional description. To obtain the renor- 
malized correlation functions in the gauge theory one must appropriately renormalize 
the bulk gravitational action using holographic renormalization. The supergravity field 
that couples to the boundary operator that we are interested in (Tr(Z 2 )) is a mixture 
of (components of the) metric on the S 5 and the four form RR potential on the S* 5 . 
It is denoted by S 22 in [28]. The one point function we want can now be obtained by 
variation of the renormalized on shell supergravity action with respect to the boundary 
value of S 22 . The result is [28] 7 



(0 S 22) =N 2 J r 3 pe 2i(p drd(f>, 



where p is the boundary condition for the LLM geometry. This boundary condition 
is most easily determined by exploiting the free fermion description of the half-BPS 
sector of M = 4 super Yang-Mills theory [2, 4]. To make the transition to the free 
fermion description, rewrite the state |$) in terms of Schur polynomials. This is easily 
accomplished using 

Ti(Z 2 ) = xm{Z) - Xfi(Z) , 
Tr(Z 4 ) = xnm(Z) - Xfin(Z) + xgj(Z) - xg(Z) , 



7 This differs by a factor of ^/2tt 2 from formula (3.60) of [28] because we are using a different 
normalization. Further, since we have dropped all spacetime dependence, we have also dropped a 
factor of e 2lt . 
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and the product rule for Schur polynomials. The Schur polynomials correspond to 
energy eigenkets of the free fermions. The energies of the free fermions in the state are 



Ei = Xi + N-i + 1 i = l,...,N, 

where Aj is the number of boxes in the ith row of the Young diagram label for the Schur 
polynomial. Using this interpretation it is straightforward to obtain the fermion phase 
space density. The density we obtain in this way is 

p = pi + cos(20)p 2 , 

where 

p 2 = —r^ — — \(Nr 2 ) N+M+2 + (N + M- l)(Nr 2 ) N+M+1 

—(N + M)(N + M + 1)(N + M - 1)(N + M - 2)(Nr 2 ) N+M - 2 

— (N + M){N + M- 1) 2 {N + M-2){N + M- 3){Nr 2 ) N+M ~ 3 ] . 

We will not need p\ in what follows. As already noted in [28], the state |<&) which is a 
superposition of a small number (for us, 6) of Schur polynomials, does not give rise to 
a smooth supergravity solution. Indeed, for a regular supergravity solution, np should 
only take the values {0, 1}; this is not the case for the above p. It is now straightforward 
to check that 

(<V0 = 



N + M 

which is the result we wanted to demonstrate. 



2.3 Multi Ring Backgrounds 

Using the results of [16, 34] we can generalize our results for the annulus to a general 
multi ring LLM geometry. The Young diagrams of these geometries are not a rectangle, 
but rather are shapes with more than 4 edges and all edges have a length of order iV 
boxes. In a geometry with m thick rings, local graviton operators were defined in 
[16]. When taking a product of the Young diagram of the background (B ring ) with the 
Young diagram of the probe (R) , one removes boxes from the probe Young diagram and 
adds them, at all possible positions, to B r i ng (respecting the usual rules for multiplying 
Young diagrams). The local graviton operators are defined by giving their product 
rule: the boxes removed from R can only be added at a specific location on B ring . This 
definition is precise enough to allow the computation of correlators; see [16] for more 
details. We can also give a rough constructive definition of the local graviton operators: 
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since we consider a boundary condition with m thick rings, the eigenvalue distribution 
in the dual matrix model will split into m well separated clumps. In the limit that 
we expect a classical geometry to emerge (large iV and large 't Hooft coupling) the 
off diagonal modes connecting these m subsectors will be very heavy and decouple. 
We expect that, when studying almost BPS states, the effect of these modes on the 
dynamics can be neglected. There is no reason to neglect off diagonal modes connecting 
eigenvalues in the same sector. Z can thus be replaced by a block diagonal matrix with 
m blocks. If clump % contains N{ eigenvalues it corresponds to an iVj x JVj block. To 
construct a local graviton operator, do not use the full matrix Z, but rather use one of 
the blocks Zj. 

It is now straight forward to verify that, if ring i has outer radius ^/N + Mj then 
we again get an exact relation 

A B {{n % ; mj }, N) = A({rn; mj}, M t + N), (2.9) 

for correlators Ab({^«; mj}, N) of gravitons localized at edge i. Thus, the gravitons at 
the edge of each ring have their own expansion parameter equal to N + M . ■ 

Finally, consider a background obtained by exciting k 2 sphere giants; each giant is 
assumed to carry a momentum k 3 < N. The corresponding LLM geometry will be a 
black disk (of area iV — fc 3 ) surrounded by a white annulus of area k 2 which is itself 
surrounded by a black annulus of area k 3 (this boundary condition is shown in figure 3 
of [23]). From the results of this section we know that excitations of the central black 
droplet could be described using either a U (N) theory with a background of k 2 sphere 
giants or using a U(N — k 3 ) theory. This is closely related to remarks made in [23], 
arrived at using a totally different approach. 

2.4 Backgrounds with > 1 Charges 

In this section we would like to study the zero coupling limit of backgrounds made 
from two matrices Z and Y. There are a number of bases which we can employ for this 
study. The basis described in [35] builds operators with definite quantum numbers for 
any desired global symmetry groups acting on the matrices. The basis of [36] uses the 
Brauer algebra to build correlators involving Z and Z' ; this basis seems to be the most 
natural for exploring brane/anti-brane systems. The basis of [37] most directly allows 
one to consider open string excitations [38, 18, 19, 20, 21] of the operator; this is the 
restricted Schur basis. These three bases do not coincide and a detailed link between 
them has been discussed in [39] . Casimirs that distinguish these bases have been given 
in [40]. Finding a spacetime interpretation of these Casimirs is a promising approach 
towards understanding the dual description of these bases. Although all three bases 
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diagonalize the two point functions in the free field theory limit, computing anomalous 
dimensions is not yet possible 8 . In particular, we do not know how to extract the 
multi-matrix BPS operators from these bases. 

In [43] a particularly simple set of restricted Schur polynomials were identified; 
these are the operators that we will focus on in this section. The operator is built 
using NM 1 Z fields and NM 2 Y fields. The restricted Schur polynomial has two labels 
Xfl,(ri,r 2 )- R * s an irreducible representation of Sn(m 1 +m 2 )- We take R to be a Young 
diagram with Mi + M 2 columns and N rows. (ri, r 2 ) is an irreducible representation of 
the SjVMi x Snm 2 subgroup. We take r\ to be a Young diagram with M\ columns and 
N rows and r 2 to be a Young diagram with M 2 columns and iV rows. We take both 
Mi and M 2 to be O(N). For this particular restricted Schur polynomial we have 

where Xri(Z) and Xr 2 (X) are standard Schur polynomials. The simplicity of this back- 
ground is a consequence of the fact that the above operator factorizes. These operators 
are not BPS or even near-BPS. Indeed, in appendix C we compute the anomalous 
dimension of these operators, which is given by 

A = N(M 1 + M 2 ) + 4XMiM 2 , A = Ng\ M . 

Thus, the results of this section can not be extrapolated to strong coupling. It is still 
interesting to ask if, at weak coupling, we can reorganize the usual expansion. 
After normalizing so that the identity has an expectation value of 1, we find 

{ W2) (XrAZ)XrAft)XrAY)XrA^)) ' 

Using exactly the same methods as were used in section 2.1, we find that the complete 
effect of the background on correlators of operators built using only Zs with operators 
built using only Z T s is to replace N — > N + Mi. In addition, we find that the complete 
effect of the background on correlators of operators built using only Ys with operators 
built using only Y^s is to replace N — > N + M 2 . Thus, these two types of observables 
admit a different reorganization of the 1/N expansion; the small parameter of the two 
expansions is not the same. 



8 For some progress see [20, 21, 41]. In particular, these articles show that at one loop, operators 
can only mix if their labels differ by the location of a single box; thus, it seems that the operator 
mixing should have a nice simple description. This is proved in the restricted Schur basis in [20, 21] 
and in the basis of [35] in [41]. 



-14- 



2.5 A Generalization of the BMN Limit for Two Charge Backgrounds 

The AdS/CFT correspondence forces us to take the limit of large 't Hooft coupling A, 
if we are to suppress curvature corrections in the string theory. It is thus a strong/weak 
coupling duality. This naive observation is evaded by the remarkable result of BMN[11]: 
the expansion parameter of the super Yang-Mills theory may be suppressed by large 
quantum numbers for "almost BPS operators" [9]. The existence of double scaling 
BMN-like limits provides a very rich class of tractable problems. In this section we 
would like to explore the proposed BMN-like sectors discovered in [8] 9 . The work [8] 
itself was an extension of [7]. The existence of these limits may be very relevant to 
understanding the gauge theory description of near-extremal charged black holes in 
AdS 5 . Using the technology we have developed, it should be possible to study these 
limits. 

We use the same two-charge background that we used in the previous subsection. 
We will however take M 1 = 0(VN) and M 2 = 0(y/N) with the 't Hooft coupling A 
large but fixed 10 . In this case, the charges J x = NM 1 ~ 0(jV§) and J 2 = NM 1 ~ 
O(iVi), whilst the difference 

A ./; J 2 1A.U..U. ~ N 

so that 11 

A — Ji — J 2 , r _i . 

77 = ~ N 2^0. 

J I + J 2 

Our background (plus open string excitations which do not modify this scaling) is a 
near BPS operator, in the language of [9]. These results are in perfect agreement with 
the supergravity arguments of [8] - the strong and weak coupling results match. The 
operator we have obtained is a "near \ BPS" state. The fact that the anomalous 
dimension is proportional to MiM 2 strongly suggests that it arises from the dynamics 
of the open strings stretching from the stack of Mi giants (described by Xn(Z)) to the 
stack of M 2 giants (described by Xr 2 (X))- 

Let us now recall a few well known facts, 't Hooft [42] has already made the beautiful 
observation that the perturbative expansion of a matrix model can be written in the 
form 

oo 

J> 2 " 2fl / 5 (A), (2.10) 

9=0 



9 We thank Shahin Sheikh- Jabbari for extremely useful correspondence, which lead to the results 
of this subsection. 

10 This is not quite the same limit as that described in [8]; in [8] the limit with gy M fixed is considered. 
n We kept the 't Hooft coupling fixed in order to obtain this scaling for i], which matches precisely 
what one obtains in the BMN limit. 
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where f g (X) is a polynomial in the 't Hooft coupling - there is no additional N depen- 
dence apart from the multiplying factor of N 2 ~ 2g . Here, plays the role of a genus 
counting parameter. In the BMN set up, one learnt that there is a new effective genus 
counting parameter g 2 = j^- replacing by studying the two point functions of BMN 
loops (each has angular momentum J). In terms of this effective genus counting pa- 
rameter, the effective 't Hooft coupling is g\ M lg 2 = QymN/J 2 - Can we repeat this 
analysis for the operators considered here? Towards this end, we set M\ = M 2 = M 
and compute the two point function 



( Xri (z)t Xr2 (Y)t Xn (z) Xr2 (y)) 



G 2 (N + M + 1) 



G 2 (N + 1)G 2 (M + 1) / 
where G 2 (n + 1) is the Barnes function defined by (T(z) is the Gamma function) 

n-l 

G 2 (z + 1) = T(z)G 2 (z) , G 2 (n + 1) = J] k\ . 

k=i 

We will now set N = aM 2 where a is a number of 0(1) in the large N limit. The 
Barnes function has the asymptotic expansion 

^G 2 (N + 1) = ^ log(AO-l log(JV)-| log(2 7 r)+C-(-l)+f; ^ l\ )N ^ 

where B 2g are the Bernoulli numbers and ('(—1) is the derivative of the Riemann zeta 
function evaluated at -1. Using this asymptotic expansion, it is clear that I 2 can be 
split into a product -Fnon-pcrt-^pcrt, where F non _ pert is a non-perturbative piece that can 
not be expanded in -j and F pert is a factor that does admit an expansion in -j 12 . This 
suggests that we should identify g 2 = ~ In the BMN case only even powers of 
g 2 appear in the genus expansion; here we can not be sure if only even powers of g 2 
appear. This depends on precisely how we factorize I 2 into -Fnon-pert-^pert- Using this 
genus expansion, we would expect an effective 't Hooft coupling 13 

A = g 2 YM M = -Lg 2 YM VN. 



Keeping this effective 't Hooft coupling fixed but arbitrarily small, one finds 

A = 2NM (l + 2A) = 2aM 3 (l + 2A) . 

12 This factorization is not unique; it will be fixed by the physics of the problem. We do not yet 
know how to do this. 

13 With this scaling, which is different to what we had above, rj does not scale with N, but can be 
made arbitrarily small. 
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This looks like a polynomial is A times some power of M. Clearly, in view of (2.10), 
A is indeed the correct 't Hooft coupling. To properly specify our state of intersecting 
giants, we need to specify the three Young diagrams which would label the restricted 
Schur polynomial of the giant system. Clearly, with our simple choice it seems that 
we have correctly obtained the correct Yang-Mills operator to describe the decoupling 
limits of [7, 8] to one loop. It would be interesting to compute higher loop corrections 
and verify that one does indeed continue to obtain a polynomial in A times M 3 . If 
extra M dependence does appear, it might be a sign that our simple guess for the 
background needs to be corrected at 0(A). This is not completely unexpected: these 
operators are not protected so one would expect corrections in general. Finally, holding 
A fixed we find the 't Hooft coupling A = gy M N = AV aN goes to infinity. Thus, even 
though we have large A the one loop correction can be made arbitrarily small. We are 
thus optimistic that, similar to what happens in the BMN limit, it will be possible to 
compare perturbative field theory results to results obtained from the dual gravitational 
description. 

The genus counting parameter g 2 — and the effective 't Hooft coupling A = 
g\ M M are exactly the parameters we would expect for the description of the large 
M limit of a U(M) gauge theory. This is rather natural: we have on the order of 
M giant gravitons whose worldvolume theory at low energy will be a U(M) gauge 
theory. This provides strong field theory evidence that the sector of the theory identified 
in [7, 8] is indeed captured by the dynamics of open string defects distributed on 
a boundstate of intersecting giant gravitons. Its natural to think that the decoupling 
limits of [7, 8] capture the decoupled low energy world volume theory of the intersecting 
giant gravitons in the same way that Maldacena's limit [1] captures the decoupled low 
energy world volume theory of N D3 branes. Our explicit computations show that the 
low energy world volume theory of the giant gravitons is weakly coupled even when the 
original J\f = 4 super Yang-Mills theory is strongly coupled. 

3. Beyond the Half-BPS Sector 

In this section we would like to determine the sector of the theory in which our reorga- 
nization of perturbation theory is valid. First, the near BPS operators we are interested 
in are BMN-like loops 

0({n}) = Tr (Y Z ni Y Z n2 Y Z ns Y ■ ■ ■ YZ nL ). 

To compute the two point correlator (0({n})0\{n})) B we start by contracting the Y 
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fields planarly. After performing the Y contractions we need to compute 



/ JjTr(Z n< Z tn< )\ 



\ i IB 



(3.1) 



We need to verify that the above amplitude admits an expansion with parameter ( - jV+M - ) 
to demonstrate that our reorganization of perturbation theory does indeed apply to the 
half-BPS and almost BPS sectors. 

We will demonstrate that a relation very similar to (2.6) holds for the leading 
contribution to the one point functions (3.1). We will then derive an exact relation. 
The first property we use is factorization (valid to leading order and if J2i n i ~ ^(1)) 



Factorization in the annulus background has been discussed in [16] and in Appendix 
A. Thus, we need only consider 



This computation can be completed in exactly the same way as the computation 
we tackled in section 2.1. We start by writing the loop of interest in terms of restricted 
Schur polynomials 



R is a Young diagram with n + 1 boxes and R! a Young diagram with n boxes. The 
expansion coefficients oir,r' are independent of N . The two point function of restricted 
Schur polynomials has been computed in [19, 37]. Using these results, we find 



The factor hooks/j/hooks^/ does not depend on N. To compute these amplitudes in 
the annulus background, we need the analog of the product rule given in Fig. 1. The 
product rule for restricted Schur polynomials [43] was used in [34] to compute precisely 
the restricted Littlewood-Richardson number needed here. Again, only a single term 
enters in the product 




(Tr(Z n Z^ n )) B = (Tr(Z n Y)Tr(Z* n Y*)) B . 




(3.2) 



(R,R') 




Xb(Z)X(R,R')(Z,Y) = gB(R,R>)(+R,+R>) X(+R,+R')(Z,Y). 



18 



+E! is obtained from +R by dropping a single box. This restricted Littlewood- 
Richardson number was computed in [34] from its definition; this involves a summation 
over restricted characters. There is however, a much simpler way to obtain this result. 
Once one has established the form xb{Z)x(r,r>)(Z,Y) = g B (R,R>)(+R,+R') X(+R,+R')( Z ,Y), 
reducing both sides with respect to Y gives 14 

c r , r >xb{Z)xr'{Z) = cr ; r,x+r>(Z) = c + R !+R/ gB(R,R>)( + R, + R>) X+R'( z ), 

where we have used the product rule of section 2.1: xb(Z)xr'(Z) = x+r'(Z) and where 
c t ,t' is the weight of the box that must be dropped from T to obtain V . Thus, 

CR,R' 

9b{r,r>){+r,+r>) — 



C +R,+R' 

This formula is exact. Although it is possible to compute things exactly, we will also 
make good use of the leading large N + M version of our results. For this reason we 
start with a large N + M analysis and then give exact results. 

Large N Analysis: To leading order in large N + M we have cr^ri = N and c + r :+ r' = 
N + M which reproduces the known result [34]. It is now straight forward to find 



C B (n,N) = {0(n)0\n)) 1 

N \ z hooks +i? f +R 



B 



Yl a R,R' \N + M) hooks+R/ f B 

(R,R') V 7 + J 

N_\ v-^ 2 hooksR f+R 
WTM ^ aR > R ' heels'- {6 - 6) 

Arguing exactly as we did in section 2.1 we find 

Cfi(n ' N) = {nTm) c(n ' N + M) ■ (3 ' 4) 

This is very similar to, but not quite the same as (2.6). Note also that (3.4) was only 
derived at leading order; (2.6) is however exact. Some insight into these relations can 
be obtained by noticing that at leading order in N we have 

(Tr (Z n )Tr (Z tn )) = N n , <Tr (Z n Z ]n )) = A n+1 . 



14 To reduce with respect to Y we take the derivative Tr ^p. A formula for the reduction of Schur 
polynomials is given in [44]; a formula for the reduction of restricted Schur polynomials is given in 
[19]. 
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Thus, as long as we restrict ourselves to the leading order, (2.6) and (3.4) can be 
restated as follows: consider an operator 0(Z, Z^) which is a product of factors of the 
form (2.2) or (3.2). Then 



(0(Z, Z ^(0^Z,^ Z y. (3,) 

At leading order, the only effect of the background is to rescale Z. This observation 
will be useful below. 

Exact Analysis: The box c RR > has a weight N + p where p is a number of 0(1). Then 

N + p hooks + ^ N + M + p hooks^ 

9b(r,r>)(+r,+r>) = N + M + p > hooks+R, = N + p hooks*, ' 

so that 

(Tr(Z n Y)Tr(Zt»Yi)) B = £ - ^ (n + Ifjf 

(R,R') * 

E hooks fi 1 f +R \ - 1 hooks fi 1 f +R 

/ , r/) hooks^ (n + l) 2 Jb ^N + M + p hooks^ (n + l) 2 / B 

\ - hooks fi 1 f +R x - hooks fi 1 

" (t^') h °° kSi? ' (n + 1)2 /b h °° ksR ' (n + 1)2 /b ' 

In the above sums, R runs over all hook Young diagrams with n + 1 boxes (diagrams 
with at most one column containing more than one box; see [34]) and R! runs over 
all possible ways of removing a box from the hook (there are usually 2 possible ways). 
The first term in this last expression is nothing but C(n, N + M). Now, it is a simple 
matter to verify that 



E hooks/; 1 _sr^ dq 1_ 1 n + 1 _ 1 
hooks# (n+ l) 2 ~~ ^ d R n + 1 ~~ n+ 1 n " n' ^ ' ' 



The first equality follows from the formula for an irreducible representation of S n : 
dx = ho " 1 | St ; the second formula follows upon inserting the explicit expression for the 
dimensions of the hook Young diagrams R and R'. Thus, the coefficient of f+R'/fB is 
1/n. Next consider 



C(n-1,N + M)= V ^I±±J^L = V ^i^. (3.7) 
(i ^ books*, n 2 /b ^ d R , n f B 
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In the above sum, Rl runs over all hook Young diagrams with n boxes and R" runs 
over all possible ways of removing a box from this hook (there are 2 possible ways). 
The coefficient of f +R > / f B in this last sum is 

^ — \ dpi" 1 1 
^ dw n n 

R" n 

This follows because we sum over all possible subductions R" of R' . This proves that 
(3.6) and (3.7) are identical and hence we find the exact relation 

C B (n, N) = C(n, N + M) - MC(n — 1,N + M) . (3.8) 

We have explicitly checked that it holds for n — 1, 2, 3, 4, 5, 6. C(n, N + M) admits an 
expansions in M l +N ; thus, we have demonstrated that the amplitudes Cs(n,N) admit 
an expansion in M x +N , albeit with the above extra M dependence in the expansion. M 
is the number of maximal giant gravitons making up the background; it is not surprising 
that there are amplitudes that have some additional dependence on M. The point is 
that all N dependence in our amplitude has, according to (3.8), been replaced by an 
N + M dependence which is perfectly consistent with what we found in the half-BPS 
sector. 

Given the above success one may ask if we can consider more general correlators and 
to work beyond the leading order. This more general analysis would tell us the class of 
operators for which our reorganization works, that is, the full class of operators that can 
be expanded in M l +N - Recall that the Schwinger- Dyson equations of a theory determine 
the correlation functions of the theory. Thus, one possible approach to our problem 
(now that the above analysis has suggested what to search for), is to demonstrate 
the replacement N — > M + iV at the level of the Schwinger- Dyson equations. This 
demonstration turns out to be straightforward. Start in the trivial background; consider 
first the simple Schwinger-Dyson equation 15 

= J [dZdtf] J- ((Z^Z^e- 8 ) . 

Performing the derivative we find 

n 

(Tr (Z n+1 Z ]n+1 )) = N (Tr {Z n Z ]n )) + < Tr (^ r ) Tr {Z n ~ r Z ]n )) . 

r=l 

15 Since we compute correlators using (2.1), to obtain Schwinger-Dyson equations that determine 
our correlators we must consider a zero dimensional matrix model with action S = Tr (ZZ^). 
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An easy computation (see Appendix A) shows that, in the background B, the above 
Schwinger-Dyson equation becomes 

n 

(Tr (Z n+1 Z j!n+1 )) B = (N + M) (Tr (Z n Z^ n )) B + < Tr ( zr ) Tr (Z n - r Z ]n )) B . 

r=l 

Clearly, the net effect of the background is to replace A — » A + M in perfect har- 
mony with (2.6). This conclusion is rather general: the net effect of the background, 
at the level of the Schwinger-Dyson equations, is the replacement A — > A + M sug- 
gesting that (2.6) should hold for all correlators of operators built using only Z and 
Z\ This conclusion is too quick. To determine the correlators, we should imagine 
solving these equations iteratively; the value of (Tr [Z n+1 Z^ n+1 )^j will thus depend on 
(Tr (Z"Z' n )) . We should start the process with the equation that follows for n = 

<Tr(ZZ t )) B = (A + M)N . 

The A multiplying (M + A) on the right hand side comes from Tr (1). This introduces 
an A dependence, which is not replaced by M + A. This result is in perfect agreement 
with (3.4). By carefully analyzing the Schwinger-Dyson equations for the half-BPS 
loops, one can verify that a similar problem does not occur for these loops, which is 
consistent with (2.6). The above departure from a pure M + A dependence is rather 
mild and one may hope that there is a simple generalization of our results. 

We are claiming that we have a new M l +N expansion parameter. Why not simply 
replace M + A — » (/i + 1) N with /j, = so that we have the usual expansion with 
\x dependent coefficients? There are (at least) two reasons to reject this proposal 

• The loop expansion parameter in the half-BPS sector is M + N ; the expansion 
coefficients have no additional M or N dependence. Since this sector includes 
gravitons, we should identify M l +N as the h for the graviton interactions. Indeed, 
a jj: expansion with /i dependent coefficients is an expansion whose coefficients 
change as N is changed, indicating additional h dependence in these coefficients. 

• Our relation (2.6) is exact and holds for any value of M. As the size of M 
changes the character of the expansion changes and it can be misleading to think 
that fluctuations are controlled by If M = 0(1), M + N can be replaced by 
N and we have the usual expansion as expected - this is the theory in the 
AdS 5 xS 5 background. In this O(jj) and the coefficients again become 
N independent. When M = O(N), M + N is itself of order A, and continues 
to control the size of fluctuations. This is the M l +N expansion we found above - 
for example the theory in the LLM annulus background. In this case, fj, = 0(1), 
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and the coefficients themselves are of order 1. When M = 0(N 2 ) we can replace 
M + N by M so that our M ^ N expansion effectively becomes a expansion. The 
correlators are dominated by contractions with the background and fluctuations 
are much smaller than in the usual expansion: they are controlled by ~ 
In this case, \x = O(N) so that the coefficients are unusually small - of size j^. 
Particularly in the last case, it would be absurd to suggest that we have a 
expansion parameter. 

This completes our demonstration that the BPS and near-BPS sectors of the theory 
admit a M l +N expansion. 

A study of the closed strings probing the LLM annulus background has been com- 
pleted in [45, 46, 16, 34]. One can describe the loops 0({n}) in terms of Cuntz oscillators 
(representing the Zs) populating a lattice (formed from the Ys) as in [47, 19, 20, 21]. 
The article [46] in particular, pointed out that the net effect of the background, on 
the one loop dilatation operator, was to rescale the Cuntz oscillators. This rescaling is 
nothing but the relation we have found in (3.5)! This rescaling could also be accounted 
for by using a new 't Hooft coupling 

Anew = g 2 YM (N + M). 

This again looks natural: we have the inverse of the effective genus counting parameter 
times g\ M - The net effect of the background on the one loop anomalous dimension 
operator (in this sector) is to replace A — > A new . At present, we are investigating the 
effect of the background on the anomalous dimension operator at two loops[48]. 

4. Discussion 

One of the results of the papers [36, 35, 37] is a new basis for the local gauge invariant 
operators of multimatrix models. This new basis diagonalizes the two point functions 
and allows an exact computation of two point (and to some extent, multipoint) correla- 
tors in the g\ M = limit. If in addition we consider near BPS correlators, corrections 
in g\ M are suppressed, so that we obtain a rather complete description of these corre- 
lators. This allows us to ask and answer a number of interesting questions, which may 
probe non-perturbative aspects of the dual quantum gravity. 

In this article we have considered precisely such a question: the dynamics of oper- 
ators whose classical dimension is 0(N 2 ). These operators are dual to states that have 
a large mass and consequently back reaction on the dual geometry is important. This 
is manifest in the fact that non-planar diagrams are no longer suppressed: although 
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they come with the usual suppressions, these are overpowered by huge combinatoric 
factors; the usual expansion breaks down. In a number of cases we have shown that 
it is possible to reorganize the expansion and have identified the small parameters that 
control these expansions. These results were neatly captured by surprisingly simple re- 
lations between correlators in the trivial background and correlators in the background 
of our heavy operator. 

Another interesting result that we have obtained, is that in the multi ring back- 
grounds and in the multi matrix backgrounds there was more than just one coupling. 
At first site this may seem puzzling, since the dilaton is a constant. To get some in- 
sight into what is going on, consider QCD. At high energies QCD is well described by 
a Lagrangian of quarks and gluons together with the coupling g\ M . At low energies, 
the coupling grows and one needs - somehow - to reorganize the theory of quarks and 
gluons into a low energy effective theory. The semi-classical objects in this low energy 
theory will be protons, neutrons, pions,... and there will be many possible coupling con- 
stants telling us how these semi-classical objects interact. The relevance of this story 
is that for the correlators we have studied, the planar approximation breaks down and 
one again needs to reorganize the theory. In this paper we have managed to explicitly 
perform this reorganization; the objects that we find that have different couplings are 
naturally interpreted as different semi-classical objects in the effective theory. 

There will be effects that are non-perturbative in the new expansion parameter. 
An example of a non-perturbative amplitude is the transition from a maximal sphere 
giant graviton (x[i N }(Z); [l N ] denotes a Young diagram with one column containing iV 
boxes) into KK gravitons with identical angular momentum J « N. The amplitude 
is given by [2, 49] 

ilog(g)-(l/( S J))log(J) 

where g — 1/N. This is non-perturbative in 1/JV. To get the corresponding amplitude 
in the annulus, using Schur technology we can again argue that we simply need to 
replace iV — > JV+M, so that the transition amplitude is non-perturbative in 1/(M+JV). 

The techniques of the article allow us to reorganize the 1/JV expansion for classes of 
observables in backgrounds described by Schur polynomials labeled by Young diagrams 
whose edges are all O(N). This does not exhaust the operators that are dual to 
classical geometries. Indeed, the Schur polynomials correspond to spacetimes whose 
LLM boundary condition preserve a rotational symmetry on the plane in which the 
LLM boundary condition is specified. This is a small subset of the possible 1/2 BPS 
LLM geometries. Further, even amongst those LLM geometries with a rotationally 
invariant boundary condition, it is likely that we have only dealt with a subset of the 
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geometries that can occur. Indeed, consider a triangular Young diagram T such that 
its longest column has N boxes, and each column has one box less than the column to 
its left. Thus our Young diagram has iV columns in total. We have tested, in a number 
of cases, that factorization of the expectation value 

_ (OXt(Z)xt(^)) 
{U)t ~ ( XT (Z) XT (Zt)) 

holds. This suggests that Xt{Z) is again dual to a classical geometry. We have not yet 
understood how to reorganize the large N expansion in this case. 

One very interesting class of excitations of an operator with a very large classical 
dimension, is obtained by attaching "open string defects" to the operator representing 
the background. Recent evidence [7, 8] suggests that these open string excited operators 
are dual to black hole micro states. Consequently, we expect that black hole microstate 
dynamics is captured in the dynamics of these operators. Our results (see in particular 
(3.5)) suggest that it may well be possible to write down spin chain descriptions for the 
open string excitations. Towards this end, we have explored the proposed BMN-type 
sectors discovered in [8] . The existence of these limits is sure to play an important role 
in understanding the gauge theory description of near-extremal charged black holes in 
AdS 5 . Taking M 1 = 0(y/N) and M 2 = 0(y/N) with the 't Hooft coupling A large but 
fixed, the charges J\ and J 2 scale as O(N^), whilst the difference 

A-Ji-J 2 i 
~ N 2 _ > o 

J l + J 2 

in perfect agreement with the supergravity arguments of [8]. Thus, strong and weak 
coupling results agree - as expected for a BMN-like limit. We have further argued that 
if we hold Mi = M 2 = 0(y/~N) fixed we find an effective genus counting parameter 
#2 = M and an effective 't Hooft coupling A = g\ M M. Keeping the effective 't 
Hooft coupling arbitrarily small, the one loop correction to the anomalous dimension 
can be made arbitrarily small, even though A is sent to infinity in the limit. This 
is similar to what happens in the BMN limit suggesting that it will be possible to 
compare perturbative field theory results to results obtained from the dual gravitational 
description. The above effective 't Hooft coupling and genus counting parameters are 
naturally identified with those of the gauge theory living on the world volume of our 
system of intersecting giant gravitons. It would be very interesting to determine the 
dimensionality of this theory: is it a 2d gauge theory as the results of [7, 8] imply? Of 
course, the operators we have considered here are presumably too simple to describe 
the black hole microstates; for that one needs to consider triangular Young diagrams 
[6]. The triangular Young diagrams have a large number of corners which implies a 
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large number of possible excitations of the operator. This is a significant increase in 
complexity as compared to the studies in this paper. 

Finally, studying the rectangular Young diagrams that we have considered here will 
allow us to construct spin chains describing BMN loops in the LLM annulus background. 
One can look for signatures of integrability for these spin chains. This provides a set 
up in which one can test if integrability survives (a class of) non-planar corrections. 
We hope to report results in this direction in the near future [48]. 
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A. Schwinger-Dyson Equations in the Annulus Background 

The Schwinger-Dyson equations provide a powerful approach to computing correlators 
in the annulus background. They are far more computationally efficient that the ap- 
proach based on cutting rules developed in [34]. The advantage of the cutting rules are 
their generality: the cutting rules work in any background. In deriving the Schwinger- 
Dyson equation, a crucial observation is that for the background we consider (recall 
that B is a Young diagram with M columns and iV rows) 



We make repeated use of this fact and consequently, the results of this appendix apply 
only to the annulus background. 

A.l Schwinger-Dyson Equations 

Start by considering 



Carrying out the derivative is straightforward, except perhaps for the term obtained 
when the derivative acts on the background. To evaluate this term, note that 



X b(Z) = det(Z) 



M 




d 



Xb(Z) 



d 



det(Z) 



M(Z"%det(Z) 



M 
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Thus, this term contributes 

M(Tr(Z n Z tn )) B 

to the Schwinger-Dyson equation. Next, focus on the term obtained by acting on the 
first Z in (Z n+1 Z^ n ) i:j which gives 

N (Ti{Z n Z^ n )) B ; 

these two terms combine to give the claimed N — > M+N replacement in the Schwinger- 
Dyson equations. Writing out all of the terms we have 

n 

(Tr {Z n+1 Z ]n+1 )) B = {N + M) (Tr {Z n Z^ n )) B + < Tr ( Z ' r ) Tl (Z n -' r Z^ n )) B . 

r=l 

A slightly more general Schwinger-Dyson equation which we found useful in the com- 
putation of correlators reads 

= J [dZdZ^}^ ((Z m Z^) tJ TT(Z n+1 ~ m )XB(Z) XB (Z^e- s ) 
which is easily seen to give 

m— 1 

(Tr (Z tn+1 Z m )Tr (Z n+1 ~ m )) B = ^ (Tr (Z^Z^-^Tr (Z r )Ti (Z n+1 ~ m )) B 

r=l 

+ {N + M) (Tr (ZtT^jTr {Z n+1 ~ m )) B + (n + 1 - m) (Tr (Z n Z^ n )) B . 
This starting point could easily be generalized to 

= J [d^^]^((Z m Zt»)y'It(Z^ 1 -™)0(Z,Zt)x B (Z)x B (^ t )e- 5 ) 

where 0(Z, Z^) is any gauge invariant operator. 

Computing correlators is now straightforward. We can obtain all correlators of 
the form ^fX^-Tr (Z 1H )Ti ((Z^) m ')^ using (2.6). Using these in the above Schwinger- 
Dyson equations, we can easily determine the correlators (J\ i Tr (Z ni Z^ '"')) which are 
of relevance for the near-BPS sector of the theory. 

A. 2 Testing Factorization 

Now that we can use the Schwinger-Dyson equations to compute correlators exactly, we 
can answer some interesting questions. One obvious question is if the annulus geometry 
provides a good background. For this to be the case, we need to have a factorization 
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of the background expectation values of gauge invariant observables. This implies that 
a single saddle point is dominating the gauge theory path integral. By the gauge 
theory/gravity correspondence, this saddle point represents a particular space-time 
geometry in gravity, that is, a classical spacetime has emerged. One nice general result 
follows from 

= J [dZdZ^^-(Z tJ [TT(ZZ^] n x B (Z)x B (^)e- s ) 

which implies that 

<[Tr (ZZX +1 ) B = (N 2 + MN + n) <[Tr (ZZ^)] n ) B . 
This recursion relation is easily solved to give 

n 

<[Tr (Z&)] n+1 ) B = l[(N 2 + MN + i) . 

i=0 

Keeping only the leading order 16 , we have 

<[Tr {ZZ ] )f +1 ) B = (N 2 + MN) n+1 = (Tr (Z Z 1 "))^ 1 , 

demonstrating factorization for these amplitudes. We can easily generalize this by 
considering 

= J [^t]^|-(Z l ,Tr(^^)[Tr(ZZt)]" XB (Z) XB (Zt) e - 5 ) 
which implies 

([Tr (ZZ t )] n+lr fr (Z P Z^ P )) B = (N 2 + MN + n+p) <[Tr (ZZ^)] n Tr (Z P Z^ P )) B . 
Once again this is easy to solve, giving 

n 

([Tr (ZZ^)] n+1 Tr {Z p Z tp )) B = f[{N 2 + MN + i +p) (Tr {Z p Z ]p )) B , 

i=0 

which becomes, at the leading order, 

([Tr (ZZ*)] n+1 Tr (Z p Z ]p )) B = (N 2 + MN) n+l (Tr (Z p Z ]p )) B 

= (Ti(Z p Z^ p )) B (Ti(ZZ^y B +1 , 

16 Of course, n and p (used below) are 0(1). 
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again demonstrating factorization. 

Above we have been careful to compute things to all orders. If we simply assume 
factorization and keep only the leading order, we get additional information about the 
leading behavior of various loops. For example, 

Tr(Z p+1 Z tp+1 ) JjTr(Z" l Z tril ) \ = (N + M) / Tr (Z p Z ]p ) JJ Tr (Z TH Z V ' H ) \ 

V i IB \ i IB 

pi 

+ yi ( Tr (^ r ^ tp ) Tr (zi n Tr ( zrHzint ) 

r=l \ i 
rij-l I 

J r=0 \ i^j 

becomes, after assuming factorization 

(Tr (Z P+1 Z^ P+1 )) B Yl <Tr (Z ni Zt n< )) B = (N + M) <Tr (Z P Z^ P )) B JJ <Tr (Z ni Z tn< )) B 

i i 

+ 5^ (Tr (Z r+p+1 Z tp Z n '- r - 1 Z tn ')) B II ( Tr (Z niZtnt )) B , 

The second term of the right hand side is subleading compared to the first term because 
(i) the first term is multiplied by (N + M) and (ii) the second term has one less trace 
in it. Thus it may be dropped to give 

(Tr {Z P+1 Z^ P+1 )) B Y[ (Tr (Z rH Z^ rH )) B = (N+M) (Tr (Z P Z^ P )) B JJ (Tr (Z n *Z^ n *)) B . 

i i 

Iterating this relation, the above result is clearly equivalent to 

(Tr (Z q Z u )) B = N{N + M) q . 

which we know is correct. 

A. 3 Testing the Cutting Rules 

In [34] a method to compute general correlators in any arbitrary LLM background was 
given. Now that we have an efficient way to compute correlators in the annulus back- 
ground we can ask: Do the cutting rule methods of [34] really work? In this appendix 
we will compute a specific correlator, first using the Schwinger-Dyson equations and 
then using the cutting rules. We have found complete agreement between the cutting 
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rule result and the result from the Schwinger-Dyson equations for any correlator we 
have computed. The cutting rules work. 
Starting from 

= J [dZdZ^^ ((Z^Z m Z^Z n+ ^ m+ % X B(Z)xB(^)e- s ) 

we obtain 

m— 1 

(Tr ^ n+1 Z m Z jfp Z n+p - m+1 )) B = < Tr {Z ]n Z r )Ti (Z ]p Z n+p - r )) B 

r=0 

n+p— m— 1 

+ < Tr (Z tn Z m Z^ p Z r )Tr ( Z n+p - m - r )) B + (N + M) <Tr ^ n Z m Z^ p Z n+P ~ m )) B . 

r=0 

Setting n — 0, m — 1 and p = 1 we have 

(Tr (Z t ZZ t Z)) B = (2JV + M) <Tr (Z f Z)) B = (2N + M)N(N + M) . 

We will summarize the cutting rule computation; for more details the reader should 
consult [34]. Evaluating this correlator using cutting rules, there are four contributions: 
the term with no contractions with the background gives 2N 3 . The terms coming from 
contracting one Z in the loop with a Z^ in Xb (Z^) give 

The terms coming from contracting both Zs in the loop with Z's in xb(Z*) give 

(tt( ±d±A\ Xb (Z) Xb (&)\ = M 2 N - N 2 M. 
\ \dZ dZ^dZ dZy AiA MiA V 

To evaluate this last contribution we had to cut the trace of four derivatives into a 
product of two traces, each containing two derivatives. This is accompanied by a 
nontrivial trace insertion factor that we evaluated in representation B. Summing these 
terms we have 

2iV 3 + AN 2 M + M 2 N - N 2 M = (2N + M)N(N + M) . 

B. Schwinger-Dyson Equations for > 1 Charge Background 

The background of interest in this appendix is (recall that r\ is a rectangular Young 
diagram with N rows and Mi columns and r 2 is a rectangular Young diagram with N 
rows and M 2 columns) 

Xn (Z) Xn (Y) = (det(Z)) M Hdet(Y)) M2 . 
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The Schwinger-Dyson equations continue to provide a powerful approach to correlator 
computations, when we consider this background built using more than one matrix. In 
this appendix we will give a few example computations. Consider the identity 17 

= J [dZdZUYdY^ ((^F m Ft mz t-i )ijXri(z)Xri(z t )Xr2(y)Xr2(y t )e -^ 
which leads to the following Schwinger-Dyson equation 

(Tr (z n y m y tm z tn )) (ri ={n + Mi) (Tr (z n - 1 y m y tm z tn - 1 )) 

n-1 

+ < Tr ( zr ) Tr (z n - 1 - r y m y tm z tn_1 )) (ri ra) . 

r=l 

If we had been working in the trivial vacuum, the only difference would have been 
to replace N + M 1 in the above equation by N. One way to think about the above 
Schwinger-Dyson equation is that to go from the left hand side to the right hand side, 
we perform one of the Z Wick contractions. To obtain the equation that follows when 
we perform a Y Wick contraction, start with the identity 

o = J [dZdtfdYdY^-jjL ((y"z m zt m yt«- 1 ) u . Xri (z) Xri (zt) Xr2 (y) Xr2 (yt) e - s ) 

which leads to the following Schwinger-Dyson equation 

(Tr (Y n Z m Z^ m Y Ul )) {ri = (N + M 2 ) (Tr {y n ~ l Z m Z^ m Y Ul - 1 )) {ri 

n-1 

+ < Tr ( yr ) Tr (Y n - 1 - r Z m Z ]m Y^ n - l )) {ri . 

r=l 

If we had been working in the trivial vacuum, the only difference would have been to 
replace N + M 2 in the above equation by N. This structure is parallel to the structure 
we found for backgrounds constructed using a single matrix: in this case we have found 
that to reproduce correlators of operators built only using Zs or Z's (Ys or y^s) we 
simply replace iV — > N + Mi (N N + M 2 ). This structure is again emerging at the 
level of the Schwinger-Dyson equations. 

Consider the last Schwinger-Dyson equation given above. In the large N limit the 
first term on the right hand side gives the leading contribution. The second term has 
one more trace in it whilst the first term is multiplied by (N + M 2 ). Naive counting of 
powers of N would suggest that these two terms are of the same order. However, the 

17 Thc action S = Tr (ZZt) + Tr (FY*). 
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leading contribution to the second term (Tr (X r ))(r l5 r 2 ) (Tr (yn-i-f^'nytn-i^ 
vanishes. Dropping the second term and iterating we find 

(Tr (Y n Z m Z^ m Y^ n )) { = (N+M 2 ) n (Tr (Z m Z^ m )) { = N (N + M 1 ) m (N + M 2 ) n . 

(B.l) 

This looks very similar to the relation (3.5); indeed, we can write 



«**. - (o v™^ 

Relations of this type would again be very useful in deriving spin chains for loops in 
this two matrix background. 

Finally, although we have already argued that there are large 't Hooft coupling 
corrections to the background, it is still interesting to ask if factorization holds. Such 
backgrounds are naturally interpreted as classical backgrounds that receive curvature 
corrections. For operators that do not mix Zs and Ys in the same trace, correlators 
factorize into a Z correlator times a Y correlator. Using the results of Appendix A 
we clearly have factorization in this case. For operators with mixed traces, a bit more 
work is needed. We will give some examples in which factorization is clear. Consider 



0= f [dZdZUYdY*]-^— ((y n+1 Z m Z^ m Y^ n ) 
< n Tr (Z^Z^) Xri (Z) Xri (Zt) X „ (Y) Xr2 (Y^)e 



which implies 

Tr (Y^ n Y n Z m Z ]m ) JjTr (Z na Z tna ) 



(r\,T2) 



= (N + M 2 ) /tt (Y* n ~ 1 Y n ~ 1 Z m Z* m ) Y[ Tr {Z Ua Z^ Ua ) \ 

\ a I ( 

n-l I 

+ Yl ( Tr ( yr ) Tr (Y n - 1 -' r Z m Z 1m Y 1n - 1 ) Y[ Tr (Z'^Z^] 



r— 1 \ a I (n,r 2 ) 

In the large N limit the first term on the right hand side gives the leading contribution 
so that we can drop the second term. Iterating, we find 

Tr (Y jn Y n Z m Z jm )Y[Tr (Z na Z^ na ) \ = (N+M 2 ) n I Tr (Z m Z^ m ) JjTr (Z na Z jria ) 

a I (ri,r 2 ) \ a 
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= (N + M 2 ) n <Tr (Z m Z^)) (rur2) J] <Tr {Z^Z^)) ^ , 

a 

where to get the last equality we used factorization of the Z, Z^ correlators. Now, use 
(B.l) to identify (N + M 2 ) n <Tr (Z m Z^ m )) as <Tr (Y* n Y n Z m Z* m )) in the last line above 
so that 

Tr (Y^ n Y n Z m Z tm ) Y[ Tr (Z na Z ]na ) 
= (Tr(Y^Y n Z m Z^)) (rur2) l[(TT(Z^Z^)) irur2) . 

a 

We can give a rather general argument for factorization: consider 

o = J [dZdzUYdY^ ((y ,H - 1 ^zt»yti')yO Xri (z)x ri (zt) Xni (y) Xni (yt) e -s) 

where O is any gauge invariant operator. This implies 

(Tr (Y^ +1 Y p+1 Z n Z^ n )0) {ri = (N + M 2 ) <Tr (Y U) Y P 'Z n Z^ n )0) 

+ < Tr (^ r )Tr (Y p - r Z m Z^ m Y^ p )0) + ( (Y p+l Z n Z^ n Y ]p ) ij -—0 ) 

r=l ri ' r2 \ ^' / (n,ra) 

The second term on the right hand side can be dropped - it vanishes at leading order. 
The third term on the right hand side can also be dropped - it represents a loop joining 
term, so that it has one less trace than the first term. If we now rescale Y — > AT N .. Y 
and Z — ■> N ^ Mi Z we recover the Schwinger- Dyson equations of the theory in the trivial 
vacuum (after dropping the same two terms justified with the same two reasons). We 
know that factorization was a property of the old Schwinger-Dyson equations so that 
we have just learnt that it is a property of the new Schwinger-Dyson equations too. 

C. Anomalous Dimension for > 1 Charge Background 

In this section we will explain how to compute the one loop anomalous dimension of 
the backgrounds considered in section 2.4. 

C.l An Identity: Excited Giant Correlators 

Correlation functions of restricted Schur polynomials have been computed in [19] (see 
also [20, 21]). The logic in these computations is first to contract the open string 
words and then to compute the remaining contractions. In this section we will obtain a 
formula that describes the result of contracting all fields except the open string words. 
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Although we derive our formula for the case of one string attached, it is simple to 
extend it to the general case. The formula we are after says 

\x { r] r {Z, W) X ( n]AZ\ W*))=A <Tr {WW*)) + B <Tr (W)Tr (W*)) . 

Recall [19] that the allowed index structure for open string word two point functions is 

(Wj{W^) = 5\5 k 3 F Q + 5}5fF 1 . 

Thus, 

( X %{Z, W) X %{Z\ = A(N 2 F + NFJ + B^F, + NF ) . 

From the technology developed in [19], we also know that 



It is now trivial to find 



/ hooks fi 2 \ f R 

= [h^ N ~ CRR,N ) N^P ' 

r ( at? , r hooks R \ f R 



C.2 Leading Contribution to Background Correlator 

The formulas we write in this section will not be general; we are considering the back- 
grounds r\ and r 2 of section 2.4. With a little extra effort one could be general. We 
are interested in computing the normalized correlation function 

\Xri {tf)Xr 2 {Y ] )XrAZ)Xr 2 {Y)) , 

to one loop. According to Appendix B of [50], the D-term, self energy and gluon 
exchange cancel at one loop order (using techniques of [51]), so to this order we only 
need to consider the contributions from the F-term. Towards this end we will now 
evaluate 

h = (XrAZ*)Xr 2 (Y*)XrAZ)Xr 2 (Y)TT([Z,Y][Z\Y*])) . 

The tricky part of this computation is the evaluation of the color combinatoric factor. 
To do this evaluation we can work in zero dimensions. Since we drop self energy 



-34- 



corrections, the F-term is normal ordered and hence the above correlator can be written 

as 



This can be rewritten, using dummy open string variables, as 



For both r\ and r 2 there is only one way possible to attach the open string. Using the 
results of the previous subsection we now obtain 

Tr ([ W ] ' w* ]) (AlTr {wwt) + BlTl {w)Tl ^ {vvi) + B2Tl {v)Tl (yt)) ' 

where 



A - Ml f A - M2 f R 
It is a simple matter to find 



h M 1 M 2 
-2NM 1 M 2 + 2- ' - 



fri fri N 

The leading contribution to this correlator comes from the terms 

Tr (ZYY^Z ] ) + Tr (YZZ^Y ] ) . 

This computation will allow us, in the next subsection, to identify and evaluate the 
leading contribution to the one loop anomalous dimension. 

C.3 Leading Contribution to the One Loop Dilatation Operator 

The 0{g YM ) contribution to the anomalous dimension 

gives A = NMi + NM 2 . To obtain the leading piece of the 0(g YM ) contribution to 
the anomalous dimension, which we have identified in the previous subsection, replace 



D 1 = -2g YM Tr : --, -- Y, Z 2g YM Tr I ZY — — J +2^ M Tr I YZ 



d d 



dY' dZ 



2 TU / TV l oJ 



d d \ „ o m / „ „ d d 



dYdZ * YM V dZdY 



The normal ordering symbols here indicate that derivatives within the normal ordering 
symbols do not act on fields inside the normal ordering symbols. It is now straightfor- 
ward to argue that 

2g YM Tr [zY^^j + 2g YM Tr (yZ-^ ^Ljj Xri (Z) Xr2 (Y) 
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= 4g 2 YM NM 1 M 2X r 1 (Z)Xri(Y). 
We are interested in two cases: 

• Both Mi and M 2 are O(N). Holding #y M A = A fixed and large (which is the 
regime in which we expect that we can trust the dual geometry), we find a one 
loop correction of 0(N 2 ) times A to the tree level value which is itself 0(A 2 ). At 
large A the tree level and one loop results are of the same order. 

• Both Mi and M 2 are 0(y/~N). Holding #y M A = A fixed and large, we find a one 
loop correction of 0(A) times A to the tree level value which is itself 0(A 3 / 2 ). 
At large A the tree level result dominates the one loop correction. 

Notice that if we take M 1 to be 0(A) and keep M 2 to be 0(1) or if we take M 2 to be 
0(A) and keep M x to be 0(1), the one loop correction becomes negligible as compared 
to the tree level value, as we would expect. It is also interesting to note that our 
operator is an eigenoperator of Di, so that at one loop and at large A it does not mix 
with other operators. 

D. Notation and Useful Results 

If R is a Young diagram, R' is a Young diagram obtained by removing one box from 
R. 

We have used the weight 18 and the hook of each box in a Young diagram. A box in 
row % and column j has a weight equal to A — % + j. We use the notation fa to denote 
the product of the weights of Young diagram R. For example 

/p = A(A + 1)(A-1). 

We use crri to denote the weight of the box that must be dropped from R to obtain 
Rl. 

To obtain the hook associated to a given box, draw a line starting from the given 
box towards the bottom of the page until you exit the Young diagram, and another line 
starting from the same box towards the right until you again exit the diagram. These 
two lines form an elbow - what we call the hook. The hook length for the given box 
is obtained by counting the number of boxes the elbow belonging to the box passes 

18 These are not the Dynkin weights. We use the terms "weight of a box" and "factor of a box" 
interchangeably. 
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through. The notation hooks^ means the product of hook lengths of Young diagram 
R. Thus, for example 

hooksrp = 1 ■ 1 • 3 . 

We have made extensive use of the exact two point function of the Schur polynomial[2] 

(xr(Z*)xs(Z)) = 5 RS f R 

where 5 R s is one if R and S are identical Young diagrams and zero otherwise. We have 
also used the product rule 

Xr(Z)xs(Z) = ^gnsTXriZ), 

T 

where g R sr is the Littlewood-Richardson number. 
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